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Abstract. Effects of vacuum polarization modify the energy levels in atoms with an orbiting particle 
heavier than an electron. The dominant effect is due to the Uehling potential. In this paper we consider 
the relativistic corrections to the energy levels caused by the Uehling potential and in particular the fine 
structure in muonic and antiprotonic atoms. We derive general expressions and consider in detail specific 
regions of parameters which allow simple asymptotic expansion. We take into account the recoil effects 
and anomalous magnetic moment in the case of an antiproton as the orbiting particle. 

PACS. 36.10.Gv Mesonic atoms and molecules, hyperonic atoms and molecules - 31.30.Jv Relativistic 
and quantum electrodynamic effects in atoms and molecules - 32.10.Fn Fine and hyperfme structure 
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1 Introduction 

Quantum electrodynamics (QED) corrections to the en- 
ergy levels of bound states have very different structure 
from one atom to another. The dominant QED effect for 
exotic atoms with an orbiting particle heavier than an 
electron is due to the so-called Uehling potential, which 
is caused by a free electron vacuum polarization loop. In 
a previous paper pQ we studied the nonrelativistic effects 
which are in particular responsible for the Lamb shift. 

In principle, the results were known for a few low-lying 
levels for a while [2]; however, they are rather complicated 
and the aim of pQ was to derive simple expressions for arbi- 
trary states, which may be realized as certain asymptotic 
expansions. For instance, the expansions are possible over 
a small parameter Za (where Z is the nuclear charge) and 
a large parameter Za m/n 2 m e (where to is the mass of an 
orbiting particle and n is the principal quantum number). 
Some other expansions are also possible. 

Analysis in pQ was performed in the leading non-re- 
lativistic approximation. Here we consider the first rela- 
tivistic correction. This correction is the dominant QED 
contribution to the fine structure and we pay special at- 
tention to the fine structure effects in muonic and antipro- 
tonic atoms. Certain nonrelativistic asymptotics were pre- 
viously studied in f4"|| 1 > ] , while in we considered a specific 
case of high n which may be of particular interest due to 
antiprotonic atoms. 

Here we derive a general expression for the dominant 
correction to the Lamb shift at a medium value of the 



nuclear charge Z . This correction is j-independent and we 
also found a relativistic correction which depends on the 
angular momentum j. Finally, the vacuum polarization 
corrections are presented in the form 

AE(nlj) = 



- (Za) 2 ^ [F nl (n n ) + (ZafH nlj (K n ) + ...], (1) 



where 



k 1 Za to 

n n Top 



(2) 
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The first term (the -F-term) is responsible for most of the 
Lamb shift and was already studied in detail in yQ. The 
relativistic correction (the ff-term) is a small correction 
to the Lamb shift, but it is responsible for the dominant 
radiative contribution to the fine structure, which we de- 
fine here as splitting for the case of j ^ j 1 while n = n' and 
1 = 1'. Due to that we consider in detail the fine-structure 
difference 

^nf( K n) = H n j J+1/2 (K n ) - H n<lil _ 1/2 (K n ) . (3) 



In the atoms of interest the mass to and parameter n n 
can be to ~ 207 to € and n n ~ 1.5 Z/n (in a muonic atom) 
and to ~ 1836 TO e and n n ~ 13.4 Z/n (in an antiprotonic 
atom) . 
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2 The relativistic expression for the Uehling 
term and the fine structure 

An exact relativistic expression for the Uehling correction 
to the energy is of the form 

AEK) = Jdrr 2 (\f nlj \ 2 + \g nlj \ 2 ) V v (r) , (4) 

where f n ij and g n ij are the upper and lower components of 
the Dirac wave function [3j and for the Uehling potential 
we use the Schwinger's parametrization 1 @ 



a f 1 , v 2 (l-v 2 /3) 

Vu{r) = - dv 

7r ./ 1 - v A 



Za 
r 



—A? 



with 



A = 



2m e 



(5) 



(0) 



VT^v 2 

The exact calculation can be performed and analytic 
results for certain states were found in [J]. Here, at the first 
stage we derive a similar expression for an arbitrary state 
in a hydrogen-like atom. Applying the known expression 
for the wave function of the Dirac-Coulomb equation (see, 
e -g-i El) an d integrating over coordinates, we obtain 



E 

i,k=0 



AE(nl ) = - 

(-\y+ k 



a 2 r(2( + n r + l)(n r )l 
— n ^ x 

n 



r(2( + i + k) 



i\(n r - i)\k\{n r - k)\ r(2Q + i + l)r(2< + k + 1) ' 



Za 



— v\ + (n r — i)(n r — k) 



E 



iiij 



Za 



v (2n r 



k) x 



where 



Kl,2,i+k+2(( K n) — g K 2 .2,i+k+2C, ( K n) 



= ( _ 1} i+m/2 (j + 1/2): 



(7) 



C = V" 2 - (Za) 2 , 



x] = ^1- (Enij/m) 2 , 

n r = n — \v\ , 

K n = nr]K n j{Za) , 

and E n ij is the exact relativistic energy of the nlj state 
for the Dirac-Coulomb problem. 
The base integrals, defined as 



K abc (K) = / dv ^W2 , 

Jo (1 ~ v ) ' \ 1 + KVl 



(8) 



1 Throughout the paper we apply the relativistic units in 
which c = h — 1 . 



are known in a closed form fjQ (cf. @j): 



_ ,'c c 1 6 c 1 3 6 c Q 

X3F2 -,- + -, 1 V -: -, o H h -; « 

1 2 2 2 2 2 2 2 2 2 



t c+1 B { a 



I 3 b c 
2' 2 ~ 2 + 2 



(9) 



c 1 3 6 c 3 be 
X3-F2 - + 1, - + -, - - - + -; -, a + 2- - + -; k 



2 2 2 2 2' 2 



2 2' 



Here 3F 2 (a, /3, 7; <5, e; z) stands for the generalized hyper- 
geometric function [7] and B (a, /3) is the beta function. 

Such expressions containing the hypergeometric func- 
tion 3F2 are cumbersome and far from being transpar- 
ent; however, various simpler asymptotic expressions are 
known ^flS^. Here we intend to expand in Za. The lead- 
ing term of the AE(nlj) expansion in Za (see the F term 
in (Tj^)) was studied by us in detail previously fjQ (see also 
[S]). Now, we derive from Eq. (QJ the leading relativistic 
correction to AE^ (nlj), related to the H term in (Tj^). 

In the limit Za — > Eq. Q transforms to an expression 
for the nonrelativistic correction to energy 



AE^(nl) = -^f(n + l)\ V 

i,k=0 

(21 + i + k + iy. 



(-iy+ k n r \ 



(21 + i + iy.(2i + k + ly. 

1 



il(n r — i)\k\(n r — k)\ 
(10) 



Ki.2,2l+i+k+2(Kn) ~ ^ -?^2,2,2i+i+A:+2 (^n ) 



This expression differs in form from those in [5] and fj^, 
but agrees with them. 



3 Results for the low lying levels 

The explicit expression for the H term is 

(n r y.(j + n + 1/2)\ (-iy+ k 



E 



2(2j + l)(n - v) ^ i\(n r - i)\k\(n r - k)\ 



(2j + i + k)\ 



(2j + i + l)\(2j + k + 1)\ 

, „n K2i + l)(4n-2j-j-fc-l) 
2(z + fe - 2i/) + ^ ^ '- 

+ (2j + 2i-2v+ l)(2j + 2k-2v + l) 

x I ip(2j + i + k + 1) + V(j + n + 3/2) 



-V(2j + i + 2)-iK2j + /c + 2)- 



(3n-2i/)(2n-2j- 1) 
4n 2 (n — v) 
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Table 1. Asymptotics of the relativistic corrections Hnljix) 
for some low-lying states. For nlj states x — K n . 



nl, 



H n lj (x) 



Is 



1/2 



i 

2X 2 



2S1/2 

2pi/2 



ZL_ _|_ 103 

Q I 70 



■eM2i) 



72 

Hi 4- 
° ' 216 



35 
16a? 

35 
Jii.i - 



3S1/2 


-§lm» 


^ 2 


+ 


25 67 
18 Tg^ 1 


3P1/2 


— |ln(2x) 




+ 


335 67 
216 18^ 


3^3/2 


-iln(2as) 


18 


+ 


8279 25 
10800 lSx' 2 



4si/ 2 
4pi/a 
4d 



4/r 



3/2 
/2 



-!M2z) 
-iln(2x) 



T34 ln (2^--. 



7T 2 _i_ 2327 _ 335 
9„ 1728 64^ 
jrf_ i 63287 _ 335 
9„ 43200 64^ 
it 2 . 16243 _ 35 
IS,," 1 " 21600 Tte 5 ' 
7rf_ , 3038863 _ 71 
27 "r" 6350400 &S 3 



x ( ^l : 2,2j"+j+fe+l(Kn) — — -fr2,2,2i+j+fc+l( K n) 

-{2j + 2i-2v + l)(2j + 2fc — + 1) 



(2n-2j - l)(2j + i + fc + 1) 



x ( -K"l,3,2j+i+fe+2(K n ) — --f^2,3,2j+i+/c+2(K„) 



(11) 



where 



dK abc (n) 
dc 



Here V"( z ) stands for the logarithmic derivative of the Eu- 
lcr Gamma function 



n— 1 

V>(n)=^(l) + £r, 



fc=i 



t/>(n + 1/2) = 2-0(2n) - ^(n) - 2 In 2 



^(l) + 21n2 + 2 £^_, 



fc=0 



and — ?/>(l) = C is the Euler constant, which finally cancels 
out in ijTTJl . 

Various high-K asymptotics of K a b c were studied pre- 
viously and applying lfTT|) to the low-lying states we derive 
results for the relativistic correction H n ij, which are pre- 
sented in Tabled for states with j = I + 1/2. Separately 
we present in Table |2] the asymptotic results related to 
the difference A^f defined in Eq. J2J. The difference is 
the leading Uehling contribution to the fine structure and 
thus it is the part of H n ij(x) of most interest. 



Table 2. Asymptotics of the Uehling contribution A^f(x) at 
x 3> 1 for the fine structure of some of the lowest states. For 
the nl states x 



til 



Alf(x) 



In - In It + - 215 , 27 



3p 



2 l n 2a; 4- si - 355 + ? 

9 ~ 18 432 T 3^3 



3d ^In2a; + Hl-^ 



4d 1 In 2a + si - 32101 , ^95 
6 111 T ig 43200 ^ 64~^ 

4d -1- In 2a; + ^ _ 36971 , _69 

in lg ill -r 54 129600 T 64~^ 

4 f i In 2a; 4- - 246373 + ^ 

3 6 111 ^ ' 108 1587600 ^ 64^ 



4 Results for the near circular states 

The general expression (|11|) contains a double summation. 
However, for circular (I = n — 1, j = I + 1/2) and near 
circular states only a few terms contribute to the sum and 
the result is relatively simple. In the limit n n 3> 1, the 
result for small values of n — I reads 



,n—l,n—l/2\l^>n 




2n 2 ,,, , 10n 3 - 21n 2 - 4n+ 12 
-ip'{2n) + - 



3(n- 1) 

2n 3 + 6n 2 — 3n + 1 n 2 
4n 2 (n-l) " 2 




(13) 



n + 3 



3(r 



,n— 2,n— 3/2(0 = 

! / \ 2n 2 

- (- ln(2« n )+V»(2n)-V(l)J - 3^3jy^(2n) 

20n 4 - 76n 3 + 7n 2 + 64n - 24 



18(n- l) 2 (2n- l) 2 



2n 3 + 6n 2 — 3n + 1 n 



1 



n + 6 



4n 2 (n — 1 ) k 2 
■ln(2«; n )+^(2n)-V'(i; 




2n 2 



(14) 



-^'(2n) 



n 2 l3(n-2)V 7 3(n-2) 

40n 6 - 316n 5 + 510n 4 + 565n 3 - 1951n 2 + 1482n - 342 
+ 18(n -l)(n- 2) 2 (2n - l) 2 (2n - 3) 
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6n — 6n + 4 n 



2n 2 (n 




The above asymptotics prove to be especially simple 

(15) 

for the case j < I: 



The natural parameter of the expansion is n/n n , not just 
1 j K n , as one can observe in a limit «>1 



(n + j + l/2)l(2j + 3)M 



In 



1 

18 



1 n 2 
2^ 



n 2 (2i + l)3(n r -l)!(2i)!(2i + 4)!!" 

(21) 

In the case of 1ow-k another asymptotics can also be 
(16) of interest, namely a result in the leading order in n but 
exact in (Za): 



This parameter can be easily understood from the presen- 
tation in the coordinate space (see PP for more discussion). 
For the fine-structure splitting we find 



,FS 



ip(2n)+ip(l)+-<ip'(2n) 




7 ? 2 J r(2C + ?v + i) 



2< 



C + i 



18(?i- l) 2 (2n- 1) 



1 n 2 

4n(n — 1) k 2 




(n r )\(^-v)r(2Q 



4C 2 (2C + 3) 
5(3/2, C)x 



Za 




(17) 



(22) 



A 



l-'S 



3 (n- l)(n-2) 



x ( ln(2 Kn ) - V(2n) + V(l) + ^ ^'( 2n ) 



6 Relativistic corrections in the logarithmic 
approximation 

To verify our calculations we consider a limit In K n 3> 1 
and find the logarithmic terms within the effective charge 



n(l6n 6 - 424n 5 + 1764n 4 - 2816n 3 + 1973n 2 - 576n + 54^ pproacn with the he i p of a substitute 



18(n - l) 2 (n - 2) 2 (2n - l) 2 (2n - 3) 



8n J - 9n 2 + 13n - 6 n 



An(n - l)(n - 2) k 2 
and its asymptotics for n»l are 

aFS = aFS 
n,n — 1 n.n- 




(18) 



Za — > Za(n n ) = Za 1 



2a 

— In « r 

37T 



(23) 



We find 



H 



In 



2U 2 



(19) 



log 
l nlj 

vFS, log 



(«n) 



1 



4 1 

3n Vi + 1/2 
4 1 1 



3 

4n 



- — In k„ , 



In «v. 



5 Asymptotic behavior at low K n 

We are mostly interested in k„ > 1 and n n 1, however, 
for the highly excited states one can arrive at the situation 
when k is large, but n n is small. 

For the case of low n n we obtain asymptotics of Eq. (|llf) 
for small K n 



3 n 1(1 + 1) 

in agreement with the direct calculations above. 



(24) 



(n + j + l/2)!(2j + 3)!! 



2(2j + l)3(n - u)(n r )\(2j + l)!(2j + 4)!! 



(2n-2j - l) 2 ^ 



+ (2j-2z/+l) 2 



In k„ + + 1/2) + + n + 3/2) 
1&(2j + 2) - + 3) - ^— + -L 



1 2j + 3 



2j - 2n + 1 
4n 2 



2j 



(20) 



7 Corrections due to the finite nuclear mass 
and anomalous magnetic moment of the 
orbiting particle 

One feature of most exotic atoms is that a gap between 
the mass of the orbiting particle m and of the nucleus M 
is not so large as in a conventional atom (where m = m e 
and thus m/M ~ 1/(1836 A), where A is the mass num- 
ber). As a result, various recoil effects become important. 
Another important property of some exotic atoms is the 
presence of the anomalous magnetic moment of the orbit- 
ing particle, which is not small. There are two kinds of 
orbiting particles with spin 1/2 which can be described 
by the Dirac equation. While the recoil effects are impor- 
tant for both muonic and antiprotonic atoms, the effects 
of the non-zero anomalous magnetic moment take place 
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only in antiprotonic atoms (k 1 = (g p — 2)/2 ~ 1.79) 2 . The 
muon anomalous magnetic moment can be treated per- 
turbatively and taken into account together with related 
radiative corrections. 

We also note that the fine structure is considered as a 
structure of levels due to the interaction of the spin of the 
orbiting particle and the orbital moment. In this consid- 
eration we neglect the nuclear magnetic moment and thus 
the hypcrfmc effects. If the nucleus has a spin, the energy 
should be averaged over it to be compared to our result. 
Such an asymmetric treatment is reasonable when the nu- 
cleus is either spinless or heavier than the orbiting particle 
and cannot be applied to, e.g., protonium, a bound system 
of a proton and antiproton. In the case of heavy nuclei it 
is customary to separate the fine and hyperfine structure 
and for this case we consider below the recoil effects in the 
fine-structure splitting. 

Below we take into account the finite nuclear mass and 
anomalous magnetic moment of the orbiting particles. To 
do that we re-visit the well-known leading contributions 
to the fine structure, found perturbatively, and show that 
their dependence on the nuclear mass and the anomalous 
magnetic moment of the orbiting particle is of a universal 
form. That allows us to adjust the results for Z\. 

Above we calculated the leading contribution to the 
fine splitting 



(25) 



in the external field approximation and neglecting the 
muon anomalous magnetic moment. 

Now, we explain how to take into account the finite 
mass of the nucleus and its anomalous magnetic moment 
for the fine structure. Similar effects are also important 
for relativistic corrections in general; however, the general 
case is much more complicated and will be considered else- 
where. 

The explicit results for A^.f(n n ) were found here by ap- 
plying the Dirac equation. Let us consider this particular 
effect perturbatively. In the case of n' = and m/M = 
we start from a Hamiltonian 



-"nr. 



where 



F7 U — 

-"nr. — 



2m R 



H FS 



+ V(r) 



(26) 



(27) 



is the complete nonrelativistic Hamiltonian for the prob- 
lem V(r) = Vc(r) + Vu(r), m R is the reduced mass, which 
in the limit m/M = is equal to the particle mass m, and 
the fine-structure effects are described by the potential 



In terms of the eigenfunctions !f{j R of Hamiltonian 



-ffwp, we present the fine structure as 



h FS = (*° K \H FS \*° 



NR./ 



(29) 



where h F s is a reduced Hamiltonian for variables related 
to spin and angular momentum only. The complete fine 
structure, which is not necessarily related to the pure 
Coulomb problem, is eventually proportional to the ma- 
trix element of 



1 



2 m 



2 ' NR I ~ — I^NR/ 



(S-I) 



(30) 



We note that if the potential V(r) depends neither on 
the mass of the orbiting particle m nor on the nuclear mass 
M, the wave function can depend only on the reduced 
mass m R because of the kinetic energy and on m e in the 
particular case of the Uehling potential. 

Thence, from dimensional analysis we can find in the 
external field approximation 3 



-(Za) 4 



Alf(K n ) = mG{m/m e ) 



(31) 



where G(x) is a dimensionless function, which depends on 
a state #$r- 

What happens if we take into account the anomalous 
magnetic moment and the finiteness of the nuclear mass? 
The nonrelativistic Hamiltonian is now of the same form 
as before (see l)27fl ) but the reduced mass 



Mm 



m R 



M - 



- m 



differs from the mass of the particle m. Meanwhile, the 
Hamiltonian for the fine structure is of the form (cf., e.g., 
[S] for the pure Coulomb case) 



H 



FS 



1 1 dV(r) 



2m 2 



di 



(8-1) 



:i + 2k') 



2m 
~M 



(1 



m 
rn 



1 



1 



dV(r) 



r m. 



Or 



(s-1). (32) 



Thus, the result reads 



mAlf( Kn ) -> -§ 
m z 



and the splitting is 



i. 2m , 
(1 + 2/0 + ^(1 



x m R Alf(K R /n) 



(33) 



H F S - 



1 1 dV(r) 



2m 2 



di 



(28) 



2 The proton anomalous magnetic moment is customarily 
denoted by n. In our paper (see, e.g., Jll9p . we use k for the 
ratio (Zam)/m e and here we denote the proton anomalous 
magnetic moment by k' . 



a 



TV 



(Za)' 



mc m 



R 



n- 



?/;■> 



(1 + 2/0 + ^(1 + *') 



3 A similar result takes place for the whole fine-structure 
splitting and for a contribution to it in any order in a once 
we consider a perturbation of the Coulomb potential by the 
Uehling potential. In our case, only a term linear in a is needed. 
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x Alf(K R /n) . (34) 

This result follows from dimentioanl analyses in (|31|l once 
we take into account all m-factors in (|32|1 . It includes the 
complete mass dependence and the anomalous magnetic 
moment and we can now use the results obtained above 
for the fine-structure difference 4 A^f (n n ). 



1(1 



+ 1) V 3tt \ n J J 



" FS c (K R \ 
2n Z A nl — 



where we introduced 



(36) 



8 Summary 

There is a wide variety of exotic atoms, which have been or 
may be successfully produced; however, only in two kinds 
of such atoms the orbiting particle has spin 1/2, namely 
in muonic and antiprotonic atoms. In this paper we study 
fine structure caused by the interaction of the spin of the 
orbiting particle and the orbital moment. Effects of the 
nuclear spin and magnetic moment are neglected. That 
is basically correct in two cases. Firstly, one can study 
an exotic atom with a spinlcss nucleus. Secondly, one can 
consider the energy averaged over the nuclear spin. The 
latter is correct as an approximation and the corrections 
enter in order (ra/M) 2 . 

Concluding, we present here expressions for the fine 
splitting in both muonic and antiprotonic atoms which 
consist of the leading term and the Uehling correction and 
properly include the anomalous magnetic moment and re- 
coil effects. 

Below we summarize results known in closed analytic 
form for antiprotonic and muonic atoms. While the kine- 
matic corrections discussed in the previous section are 
well-known, in this paper we found various presentations 
for the leading radiative correction (A^f ) in Sects. 3-5. 



8.1 Antiprotonic atoms 

The QED result for the fine splitting in antiprotonic atoms 
is of the form 

at? i n 1 {Za) 4 m R c 2 /m R \* 
AE FS (nl) = 3 I 

2m , 

(1 + 2k') + — (1 + k') 



Alf' c (x)=A-f(x) 



PS / 



1 



1(1 + 1) 



\ n ) 



(35) 



or 



at? i n 1 (Za) 4 m R c 2 fm R y 

AE FS (nl) = - g ( 

2 n 6 V m / 



x 



/n 2m , 
(1 + 2*0 + ^(1 



K 



4 We apply here the dimensional analysis only to a term 
linear in perturbation Vu', however, the scaling factor, which 
takes into account the nuclear mass and the anomalous mag- 
netic moment of the orbiting particle, is correct in any order 
and for any nonrelativistic potential, as long as the latter does 
not explicitly depend on m or M. 



A FS, log 



which does not contain the leading logarithmic term (|24|l 
for high n n . The higher order corrections include a small 
factor which is either a or (Za) 2 . 

The QED result above does not include any nucleus- 
related effects such as the finite-size effects, virtual QED 
annihilation effects and effects of strong interactions. The 
latter may be dominant for low-^ states. On the contrary, 
for higher-Z states the dominant correction to the fine 
structure is determined by the equations above. 



8.2 Muonic atoms 

The QED result is of the form 



„ ^ , ,n 1 (Za) 4 m R c 2 
AE FS (nl) = - K - 



7T \ m J 



1 



7T V n I 



AE FS (nl) 



\ M J 7T V m / 



1(1 + 1 

1 (Za) 4 m R c 2 r 
2 

W+-T) ( 1 + ^ ln (v) 



(37) 



(38) 



The higher order corrections include a small factor 
which is either a or (Za) 2 . Some of them have been known 
only for a few particular cases, but not in general. The 
higher order effects for some particular situations are re- 
viewed, e.g., in |1(JII11) . 

The difference with antiprotonic atoms is a perturba- 
tive treatment of the anomalous magnetic moment of the 
orbiting particle. The accuracy of the muonic expression 
is higher than in the antiprotonic case, because of lack 
of the direct strong interaction and smaller importance of 
the nuclear finite-size effects because of a larger average 
distance between the orbiting particle and the nucleus. 
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